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Abstract 

Gravitational corrections in J\f — 1 and J\f — 2 supersymmetric gauge theories 
are obtained from topological string amplitudes. We show how they are recov¬ 
ered in matrix model computations. This provides a test of the proposal by 
Dijkgraaf and Vafa beyond the planar limit. Both, matrix model and topolog¬ 
ical string theory, are used to check a conjecture of Nekrasov concerning these 
gravitational couplings in Seiberg-Witten theory. Our analysis is performed for 
those gauge theories which are related to the cubic matrix model, i.e. pure 
SU( 2) Seiberg-Witten theory and Af — 2 U(N) SYM broken to Af — 1 via a 
cubic superpotential. We outline the computation of the topological amplitudes 
for the local Calabi-Yau manifolds which are relevant for these two cases. 


November 2002 


1. Introduction 


The fact that string theory can be a powerful tool to study four-dimensional supersym¬ 
metric gauge theories has been long appreciated. One of the most successful approaches 
is geometric engineering, where the gauge theory is appropriately embedded in a local 
Calabi-Yau compactification of type II string theory. A large and interesting class of su¬ 
persymmetric field theories can be geometrically engineered in this way, including Af — 2 
supersymmetric Yang-Mills theory with or without matter [I]][2j]. The gauge theory is re¬ 
covered in a certain singular limit of the Calabi-Yau geometry where the string corrections 
decouple, and the gauge group and matter content of the field theory depend only on local 
properties of the Calabi-Yau space near the singularity. 

Type II string theory on Calabi-Yau geometries has a topological subsector called 
topological string theory. It turns out that amplitudes in topological string theory com¬ 
pute exactly certain holomorphic terms in the effective action |>|. There are two types of 
topological strings, the A-model and the B-model. In the A-model, the relevant correlation 
functions are sums of world-sheet instanton contributions of a given genus. Using mirror 
symmetry one can relate the A-model to the B-model, and the latter provides an elegant 
and efficient tool to perform these sums exactly. From a pure field theory point of view 
the exact holomorphic information comes from summing over space-time instantons. The 
decoupling limit relates string world-sheet instantons to field theory space-time instan¬ 
tons, and for example the space-time instanton corrected exact gauge coupling of Af — 2 
Seiberg-Witten theory is recovered in that limit from the genus zero topological string 
amplitude. Higher genus world-sheet instanton corrections are believed to calculate ex¬ 
actly, in the same limit, certain (holomorphic) couplings of the gauge theory to gravity. 
On the field theory side, these couplings haven’t been much studied. For Af — 2 theories, 
the genus one coupling can be extracted from the low-energy twisted theory, where the 
coupling to gravity plays an important role [|6j,jj],[|||] . More recently there has been impor¬ 
tant progress in calculating also the space-time instantons contributions for higher genus 
0 . 0 . 

In the context of geometric engineering of gauge theories, one can break Af — 2 to 
Af — 1 by adding space-time filling D-branes wrapped on internal supersymmetric cycles, 
or in a dual picture by turning on RR-fluxes. The best studied holomorphic Af — 1 quantity 
is the superpotential, which is computed via genus zero open or closed topological string 
amplitudes, respectively. As in the Af — 2 case, higher genus amplitudes give gravitational 
corrections to the gauge theory effective action lfT|.[T2|| . 


I 


Recently Dijkgraaf and Vafa 


15] have made the exciting proposal that these 


holomorphic terms in Af — 1 and Af — 2 supersymmetric field theories are calculated by 
matrix models. The genus expansion is replaced by an expansion in 1/N, N being the 
dimension of the matrices. This leads to the possibility to use matrix models techniques to 
study supersymmetric gauge theories as well as topological string theories. In particular in 
the Af = 1 context the tree level superpotential appears in the action of the matrix model 
and the exact effective low-energy superpotential is recovered in the planar limit. This has 
been checked explicitly in various examples, see for example [ TS,TS[ . Although the original 
proposal of Dijkgraaf and Vafa applies to Af = 1 models, one can recover results for the 
Af = 2 theory in an appropriate limit OEM • 


The aim of this paper is to bring the above ideas together in order to perform a non¬ 
trivial check of the matrix model proposal beyond the planar limit, calculating thereby the 
corresponding gravitational terms in the supersymmetric field theories with completely 
independent methods. In this paper we will consider the simplest case of the cubic matrix 
potential. According to Dijkgraaf and Vafa, this matrix model captures the holomorphic 
couplings of IIB theory on a Calabi-Yau geometry with fluxes that generate a superpo¬ 
tential. On the other hand, Af — 2 SU (2) Seiberg-Witten theory is recovered at a special 
co-dimension one sub-locus in the moduli space, where the Af — 2 supersymmetry is re¬ 
stored. This in turn can be geometrically engineered by decoupling the string correction 
from the dual type IIA string on the canonical line bundle over P 1 x P 1 = Fo [jJJ. The 
interrelations between these theories are shown in the diagram. 


N=2 local CY string geometry 
0(-2 -2) — F 0 
world-sheet instantons 

a’—o 

V 


N=2 SU(2) SYM coupled to SUGRA 
space-time instantons 



N= 1 local CY geometry with fluxes 
uv=W / (x)“+ f (x) + y 2 


D-V 


cubic matrix model 

wt m 2 g 3 
W= — x +— x 
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We are able to check non-planar contributions of the cubic matrix model against the 
holomorphic gravitational couplings in various ways. First, we calculate the exact genus 
one free-energy F ( 1; of the matrix model over the whole J\f = 1 moduli space. This matches 
the genus one contribution in the corresponding IIB model, which we obtained using the 
holomorphic anomaly equation. Second, we consider the holomorphic gravitational cou¬ 
plings of the J\f — 2 theory. We compute them in the string decoupling limit of type IIA 
theory on local Fq. and we check these in this limit against the corresponding space-time 
instanton calculations. Finally, we show that, by specializing the matrix model answer for 
F l 11 to the J\f — 2 subspace, one recovers the right answer. 


The paper is organized as follows. In section 2, we analyze topological string theory on 
the relevant local Calabi-Yau geometries: local Fq, which geometrically engineers J\f — 2, 
SU( 2) Yang-Mills theory, and the geometries that engineer J\f = 2 U(N ) super Yang-Mills 
theory broken down to J\f = 1 by a tree level cubic superpotential W($). In section 3, we 
derive an expression for the gravitational correction to F ^ of Al — 2, SU( 2) Yang-Mills 
theory from topological held theory, and we show that it agrees with the decoupling limit of 
the corresponding amplitude on local Fo. The geometries considered in section 2 are both 
captured by the cubic matrix model, which we analyze in section 4 from a perturbative 
point of view and also from the point of view of the loop equations. In particular, we derive 
an exact expression for F { 1 ■' that is shown to agree with the one computed in section 2 for 
the Af — l theory. We also analyze the embedding of Af — 2 in the matrix model and we 
reproduce the gravitational coupling of SU (2) Seiberg-Witten theory from the genus one 
free energy of the matrix model. In section 5 we focus again on SU{ 2), Af — 2 and we 
show that the results obtained in this paper agree with the direct instanton computations 


of [10,11], confirming in this way a conjecture of Nekrasov. 


2. Topological string in the associated local Calabi-Yau geometries 


In this section we will consider topological string theory on two different pairs of non¬ 
compact Calabi-Yau three-folds I and IF They are both related to the cubic matrix model 
in a way to be described below. 

The first pair of Calabi-Yau three-folds, denoted by /, is a mirror pair. The A-nrodel 
geometry of / is the total space of the anti-canonical line bundle 0{— 2, —2) —> P 1 xP 1 . 
The B-nrodel is its local mirror geometry, which is (in a patch) given by the constraint 


£3 t 1 £3 ^2 

vw = 1 + Yi + —— +Y 2 + —— 
*1 Y 2 


( 2 . 1 ) 
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where v , w E C and Y \, Y 2 E C*. t t are the two complex structure deformations. A 
canonical parameterization of the complex structure moduli space is in terms of periods f j 
of a meromorphic differential on the elliptic curve obtained from m by setting vw — 0 [jlj]. 
The relation between tj_ and t t will be given below. The Picard-Fuchs equations as well as 
explicit expressions for the periods in terms of various expansion parameters can be found 


m 


The Calabi-Yau three-folds in the second pair, which will be denoted by II and II, 
are related by a geometric transition. The II geometry is obtained by deforming the 
bundle 0(— 2) © 0(0) —■> P 1 , such that the location of the P 1 section in the 0(0) direction, 
parameterized by x, is restricted to the two critical points, a\,a 2 , of a cubic potential 
W(x) = yX 2 + §x 3 jnj- Consider now type IIB string theory on this Calabi-Yau three¬ 
fold. As pointed out in j2(||l9||, by wrapping N\, N 2 D5 branes around the Pi’s located 


at ai, « 2 , respectively, one can geometrically engineer an J\f — 1 theory. This theory is 
U(N) J\f — 2 Yang-Mills theory broken to J\f — 1 via the addition of the superpotential 
W(5>) for the Af — 1 adjoint chiral superfield $ that is part of the Af = 2 gauge multiplet. 
The configuration with the above distribution of branes corresponds to a classical vacuum 
where U(N) — > U(Ni) x U(N 2 ). At low energies the SU(Ni) C U(Nf parts confine and 
the unbroken gauge group is a product of U( 1) factors. For each SU(N t ) one has, at 
low energy, a glue-ball superfield Si. The resulting dynamics is governed by the effective 
superpotential for the glue-balls, lF e n'(F,;). 


The II geometry is obtained from II by a geometric transition ||TT|, in which the two 
P 1 \s are contracted to a point and the resulting singular space is smoothened by two .S' 3 ’s 
of finite size. This has the local description as a hypersurface in C 4 : 


vw — W'(x) 2 + fi(x) + y 2 , 


( 2 . 2 ) 


where x,y,v,w E C and /1 is a polynomial of order one, which splits the double zeros 
of W'(x) 2 to af,ci±. The two complex structure deformations of (|2.2|) are the two pa¬ 
rameters of /1 or, alternatively, the differences af — a~. The periods S; = -fp, f a ~ w and 

Ffj = f a + * = 1,2, where u> — d x^/W'(x) 2 + fi(x), emerge by integrating over two 

i 

dimensions of the period integrals of the holomorphic three-form of the local Calabi-Yau 
geometry ( 272) ]T^] . They are functions of the complex structure parameters and of the 
parameters m and g which appear in W. The Ilj also depend on a cut-off A which must be 
introduced since we are working on a non-compact curve (i.e. x E C rather than x E P 1 ; 
otherwise there would be e.g. a linear relation Si + S 2 = 0). After the geometric transition 
the D5 branes disappear, and we are left with a closed IIB string geometry with fluxes. 
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2.1. Topological amplitudes on 0(— 2, —2) —> P 1 x P 1 and their field theory interpretation. 

The geometry / was considered in [I] to geometrically engineer SU(2) Seiberg-Witten 
field theory from string theory. The two moduli of the A-model are the complexified 
Kahler parameters t [ and 0 of the two LP 1 's. String corrections to the gauge theory can be 
decoupled by sending e ~ yfod —> 0. This should be done in such a way as to preserve the 
renormalization group relation t\ 1/fls ~ l/fi'YM rs./ log(— ^-) where mw ~ a ~ t2, A is 
the scale of the gauge theory, g s the string coupling and a parameterizes the expectation 
value of the adjoint scalar field in the Cartan subalgebra of SU(2). This leads to the double 
scaling limit 0 e — > 0 [pT|,|lH 


exp ( —\ ) = exp(—ti) = cie 4 A 4 
9s 


(2.3) 


t2 — ec2a , 


If we choose c\ = \ and c 2 = 2, A and a turn out to be those of the Seiberg-Witten theory. 
We will make this choice in the following. 

Our aim is to solve the topological /i-model on this geometry at higher genus and 
to take the field theory limit ( |2.3| ) of these contributions. In sect. 3 we will derive the 
genus one result directly as the holomorphic coupling of the J\f — 2 chiral multiplet to 
gravity and in sect. 4 we will recover it from non-planar contributions to the free energy of 
the cubic matrix model, evaluated at the minimum of the effective superpotential. It has 
been shown explicitly in |l(|] that the effective superpotential is computed by the planar 
diagrams. Our results thus provide a test of the Dijkgraaf-Vafa conjecture beyond the 
planar approximation. 

The leading structure of the period vector in the large complex structure variables 


near zs = e ti — 


is 


n = (Mi = log(>i) + a, t 2 = log(2 2 ) + cr, F q = tpt 2 + pY 


(2.4) 


where a — 2(zi + Z 2 ) + 3 (zf + ^Z\Z 2 + z |) + 0(z 3 ) and p — A(z\ + z-Y) + (9 z\ + 32z±Z2 + 
9 z$) + 0(z 3 ). The formulae in [[22[ give the whole expansions. We denote q t — exp (—ti). 
With q — qi/q 2 and Q — (p 2 we can obtain the genus zero prepotential by integrating 
F q — —2q-^F(°\q, Q ), cf. |2^]. For later use we define Q =: e~ T and q =: 


: e 


—t 


1 In the heterotic-type II duality t\ is identified with the size of the base of the K3 fibration 

and t 2 with the size of a 2-cycle in the fiber. 


A schematic view of the moduli space can be found in [18]. 
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The topological amplitudes for g > 0 are not holomorphic in the complex struc¬ 
ture parameters, but various holomorphic limits exist. For example, if we expand them 
around the large complex structure point where Qi —> 0 one can take the holomorphic limit 
qi ^ 0 ||23|| in which F^ has the form 


OO 

= c« + (2.6) 

m,n>0 

0,0 


where C <9) stands for classical terms. The expansion parameters Cnm are the Gromov- 
Witten invariants for the maps of genus g world-sheets with bidegree (m, n), while the 
classical terms arise from maps of bidegree (0, 0). In our case we have = ^t\ — \t\t 2 — 
+ Mi; C (1) = ~ It and c(9) = (~ 1 ) g 4 g (fg-f) 2 (2g-i)! X (x is the ‘regularized’ 

Euler number of the non-compact target space). 

In the decoupling limit ( |2.3|) the topological amplitudes F (jl> compute the holomorphic 
part of the low-energy effective J\f — 2 SUSY field theory. In particular the complexified 
effective gauge coupling (r = -p? + is obtained in the field theory limit as r(u) — 
— li me ^ 0 4 d^F^F [jl|], where the log(e) terms, which appear in this limit, have 
been absorbed in the bare coupling and u(a) parameterizes the field theory moduli space. 
Similarly, the holomorphic functions F^iu^a)) , which multiply the following combinations 
of the self-dual part of the graviphoton field strength, F + , and of the curvature tensor, 
_R_l_, in the effective Lagrangian 


F {9 \u(a))Fl 9 ~ 2 R 2 + 


( 2 . 6 ) 


can be calculated from F'' a! in the limit ( |2..'^ ). On dimensional grounds, these terms are 
suppressed by M^ _2p) ~ e ^ 9 ~ 2 \ 


A non-vanishing n-space-time-instanton contribution to JFO) (u) = lim^o e ( 2 9 _ 2 ) F (9) 
arises from the infinite number of terms in F^ at fixed degree n (i.e. w.r.t. the base) 
in (|2.5|) . This contribution must have the structure a 2 ~ 2g • This requires that the 

corresponding term in /ho scales like ^ ae y g - 2 ( j ~ ( 1 _ (/ 2 ) 4 n+ 2 g -2 - This in turn implies 
that Cnm grows with m as 


da) — 


= 7 nm 


4n— 3+2g 


(2.7) 


A similar argument was given in jlj for g — 0. The proportionality factor j n is directly 
related to the gauge theory space-time instantons. E.g. for the gauge coupling we have 


fi(°) 

^nm 


2(4n —1) 

Fa - FT ^ 171 

(4 n — 3 )!n 


(2.8) 
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with T r the n-instanton contribution to as calculated in [E5 . 


The topological amplitudes for this model in the large volume limit have been cal¬ 
culated for g — 0 in [|l| and for g = 1 in & Generally, the amplitudes F' !l> for g > 1 
are defined recursively in terms of all amplitudes with lower g and the propagators S' :I 
as defined in ||. It can be shown that for this model the propagators can be chosen to 
vanish except for S u |[1T| . It can be determined from the g — 0, 1 amplitudes by using the 
simplification of the holomorphic anomaly equation which occurs at genus one in the case 
with only one non-zero propagator: 


d t F^ = S u d?F^ + d t J2a r log(A r ) . (2.9) 

r 

Here the sum is over the components of the discriminant loci which are (i) are the conifold 
locus Ai = 1 — 8(21 + Z 2 ) + 16(^i + z$) — 32^1^2 and (ii) the divisors A 2 = Z 1 Z 2 at large 
radius limit. The coefficients a r parameterize the holomorphic ambiguity. A convenient 
choice is — Tjh and a 2 — which yields, using the genus 0, 1 results, the expansion of 
S u 


S U = ^ + 4qiQ 2 + 16(q?q 2 + QiqI) + 40 (qfq 2 + qiql) + 188 q\ 


'q\ + 0(q l 


( 2 . 10 ) 


With the recursive definition of the F ^ [^f, worked out for the local H-model in ^5[, the 
higher F ifl ^ can be calculated up to the holomorphic ambiguity. We have fixed the latter 
from the knowledge of the absence of holomorphic curves of low degreep 6 . 


We have calculated the Gromov-Witten invariants using the above procedure and 
found agreement with genus 2 results of [^7|, who used the same method and at higher 
genus with the ones evaluated in |2f|] using Chern-Simons theory. In particular we find 
that the cifm are 


2n —1 

0 ( „t = -Pi 4 "- 3+2s V) = p‘ 2 "" 2 + 2 »>(x) n (* + m) (2.11) 

k= 1 

where P^\x) , pif\x) are polynomials of degree d. The remarkable fact that the Gromov- 
Witten invariants lie on polynomials is a consequence of the embedding of the SU( 2) 
J\f — 2 supersymmetric gauge theory in the string geometry and holds therefore also for 
other geometries like O(-K) —> F n with n — 1 , 2 , which allow for such an embedding [[Ij]. 
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We can give the topological string amplitudes exactly to all orders in q 2 by specifying the 
Pn 2n 2+29 \x ) or alternatively by writing 


OO yi 

F«» = C (0> (t) - 2(U 3 („ 2 ) + £ -^L— h m {q ^ 


n= 1 

1 OO ryi 

1 Qi 


F< 1 )=C’W(<)--(Li 1 fe)-^ 


“ "1 1 - ®) 


4 n 


^n } (92 


( 2 . 12 ) 


_p(s) — (7(9) _|_ ( l) 3 -®2 g 


E 


<?r 


g(2g - 2)! ^ n 3 “ 2 s(l - Q 2 ) 4ri+2s_ 


32^(9 2 


The hii 11 are polynomials of degree 2(n + <7 —1) whose symmetric coefficients determine 
all orders in g 2 for given a order in gi. E.g. for g — 0 we find 

*<°> =1 

—1 + 18 q + q 2 

hf ] =1 + 98 q + 450q 2 + 98g 3 + q 4 

4 0) =1 + 306 q + 4851 g 2 + 13188 q 3 + 4851 g 4 + ... (2.13) 

/4 0) =1 + 732 q + 26903 q 2 + 206434 q 3 + 426060 q 4 + ... 

4 0) =1 + 1490 q + 105315 q 2 + 1660604 q 3 + 8358292 q 4 + 14651604 g 5 + . .. 


The missing terms, indicated by ..can be recovered from the ones displayed via 
hn\l/q) — g-( 2 ( n +9- 1 )/ l (9)(g). ( p. 13|) and ( j2T3p allow to determine the T n and also pro¬ 
vides an expansion scheme for the string corrections. E.g. using c \ = \ and c 2 = 2 we ex¬ 
pand in the limit (PT5|), using Li 3 (l—ea) = £(3) — z ^ae+(| — — \ log(ea))a 2 e 2 +(!}(e 3 ). 

This yields to the relevant order, namely e 2 , from ( |2.12|) and ( |2.13|) 


,/ . . x 5x 2 3x 3 1469s 4 4471s 5 40397s 6 441325s 7 ^ , 

r<°> = ^(logW+cl 0 *-^-^-^-33,-5TW-335-_^ +0 (V 


2 18 2 31 


5 • 2 34 


2 43 


7 ■ 2 47 


(2.14) 

where s = (-y) 4 and c^ 0> is related to the bare gauge coupling into which a log(e) term has 


been absorbed. This matches the Seiberg-Witten prepotential |24 














The corresponding expressions for genus one are 


h? = (1 - qf 

h { 2 ] = 1-2 q- 94 q 2 - 2 q 3 + q A 
h { p = 1 — 1137 q 2 — 3872 q 3 — ... 

= 1 + 4 q — 6818 Q 2 - 72168 q 3 - 158262 q 4 - .. . 

/4 1} = 1 + 10 q- 28440 g 2 - 643440 q 3 - 3622665 g 4 - 6479092 q 5 - ... 

/£> = 1 + 18 q - 94008 q 2 - 3827252 q 3 - 41834673 q 4 - 167100606 q 5 - 265697392 g 6 - ... 


(2.15) 

The holomorphic quantities of the gauge theory from higher J r(r,} ( |2.6|) can be extracted 
in the limit (|2.3|) at order e 2_2ff , as explained above. For T { 11 we obtain from ([2. 1 2|) and 


(EH) 


i 


^ (1, =-^ lo gW+C (1) + ^3 + 3 . 2l4 


X " 


1647x 4 981x 5 450137x 6 45111x 7 ^ 

+ + + ^ ^ +0(x 8 ) , 


2 29 


2 31 


3 ■ 2 41 


2 42 


(2.16) 

this expression matches exactly the one calculated in the twisted J\f = 2 gauge theory in 
the next section and the matrix model calculation in section 4. 

We have calculated the h\ 9 ' with the topological B-model methods described in fl8|1 , 
albeit at another expansion point, up to genus 3. These expressions can be found in 
appendix A. In the gauge theory limit we get the following expressions for the holomorphic 
functions in ( |2.6|) 


+”( 2 ) 

T {3) 


1 / 1 llx 2 117x 3 171201x 4 1919923x 5 

a 2 V 480 ~~ 2^ 2^ 5 • 2 37 

1/1 7x 2 293x 3 985823x 4 4069345x 5 

8064 + 2 19 + H 23_ + 2^ + 5 • 2 38 


+ 


96877135x 6 

^47 

416333277x 

2^ 


+ 0(i 7 )) . 

- + 0(x + . 

(2.17) 


2.2. The B-model on the local geometry II 


As explained before, the geometry II describes the low-energy dynamics of U(N ) 
super Yang-Mills theory with a cubic tree-level superpotential, and Dijkgraaf and Vafa 
conjectured in [13|] that the amplitudes of the topological B-model on this geometry, which 
correspond to the holomorphic couplings of the gauge theory, is captured by the cubic 
matrix model. At genus zero this was checked by noticing that the planar solution of the 
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matrix model reproduces the special geometry of II [13], and it was also confirmed in [16 
by using matrix model perturbation theory. To check this at genus one, we shall calculate 
As explained in |23[,|3|1, F n > is a section of a determinant line bundle over the complex 
moduli space, which was studied by Ray and Singer. F ^ satisfies a holomorphic anomaly 
equation which can be readily integrated up to a holomorphic ambiguity to be discussed 
below. 

To proceed we need the periods Si, S 2 from 0 . whose notation we adopt. It is 
convenient to change variables (a]~, , a^~, aj) = (x 4 , x 2 , x 3 , x 4 ) —> (A 2 i, A 43 , Q , I) where 


A 21 = ^(x 2 - xi) , A 43 = ^(x 4 - x 3 ) 

„ 1 , . m 

Q = -(x 4 + x 2 + x 3 + x 4 ) =- 

2 g 

1 = ^[(^3 + x A ) - (x 4 + x 2 )] = \/A 2 - 2(A| X + A 2 3 ) 


(2.18) 


We will also use z\ — A 2 3 and Z 2 — A^ 4 and A = («i — a 2 ) = y. 

The periods Si(zi, z 2 , g, A) were computed in H]i. To express the R-model ampli¬ 


tudes in terms of the Si we need the inverse relations 


(2 Sf - 3S 1 S 2 ) + —r V S 1 (5S 1 - 135 2 )(45 4 - 3 S 2 ) + 0(S 4 ) 

g 6 A 7 


z i{Si, S 2 , g, A) ^a^ 1 ^”^ 2 A 4 
z 2 (S 1 ,S 2 ,g,A) = z 1 (S 2 ,S 1 ,g,-A) 

(2.19) 

According to [Q and 0 and taking the simplification in the local case [[TJ into account, 
we expect the holomorphic Si —> 0 limit of the genus one B-model amplitude to be, up to 
an additive constant, 

F^ — i log (det 


dzi 


dS„ 


m 


( 2 . 20 ) 


For the holomorphic ambiguity f(z) we make the Ansatz f(z) — (J^ A^)/^ 3 . The Aj 
are components of the discriminant of the Riemann surface. It splits into two factors, the 
conifold divisor A 4 = z\-z 2 and A 2 = (x 3 — x 4 )(x 3 — x 2 )(x 4 — x 2 )(x 4 — x 4 ). The exponents 
Ki determine the leading behavior of F ir> at the degeneration loci. At the conifold we 


expect |2fJ K\ — — A We fixed the by comparing with any three coefficients of the 


matrix model computation, which we will present in sect. 4, to K\ = — |, k 2 = § and 
Ks — —1. Note that I — 0 is not a discriminant component of the local geometry. The 


3 Note that the sums in the expressions for the periods Si and IR in appendix B of [19] should 
read £~ , c„ Aand £“,i A^ 


‘ 2nG _ 2 ffl , respectively. 


i21 (2n—2)! 
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above ansatz for / was motivated by simplicity. I seems to appear in / as a branch locus 
in the moduli space similarly as the orbifold point 4/ = 0 appears in the ambiguity of 
F W for the quintic in 0. F W has the expansion 


F ^ =-log 

12 & 


g 6 

+ ^9 


m 


+ 


m 


12 


+ 


9 


10 


m 


15 


+ 


9 


12 


m 


18 


S\S 2 

A 6 


332 

~9~ 


S 


+ |(5l - s 2 ) + ^ ( -S'f - ^5 ± S 2 + -5. 
923 


m L 


sjs 2 +... 


g 8 (1864 


l 


3 

/54416 




47083 

6 


SfS 2 + 15349Si 2 S 2 2 - ... 




Sf - 187528 Si 4 S 2 + 570066Si 3 S 2 2 


/1762048 
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Sf 


12980560Si c5 c 54863776 4 2 
- —-£qh 2 H-—-<->i<->2 


256344964 


S?S 8 + ... 


+ 0(S 7 ) : 

where .. . means antisymmetric completion. 


( 2 . 21 ) 


3. F (l> from topological field theory 


In the decoupling limit ( |2.3|) , the amplitudes give corrections to the low energy 
effective action of Af — 2 super Yang-Mills theory of the form Q2.6|) . In particular, T { 11 
gives a term 


d A xF {1 \a)FvBz 


which can be integrated to obtain 




(3,1) 


(3.2) 


Here a is a coordinate in the moduli space of Af — 2 super Yang-Mills, and x and a denote 
the Euler characteristic and signature of the four-manifold, respectively. A coupling to 
gravity like ( |3.2|) is natural when the theory is embedded in string theory, as we are doing 
here. Another way in which couplings to gravity are relevant is when the Af — 2 theory 
is topologically twisted. Indeed, as shown in 0, the low-energy effective action of twisted 
Af — 2 super Yang-Mills on an arbitrary four-manifold contains terms proportional to the 
Euler characteristic and the signature (these terms appear in any twisted gauge theory with 
Af — 2 or Af — 4 [0 supersymmetry). As pointed out in [jf|,[l4| , when the manifold is 
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hyperKahler, the twisted and the physical theory agree, and we should therefore be able to 
compare the coupling ( |3.2| ) with the corresponding result in topological field theory. This 
was done in |L4]] for the J\f — 4 theory. 

The geometry I of the preceding section engineers J\f — 2 Yang-Mills theory with 
gauge group SU( 2) in the limit (|2TTJ) . On the other hand, the coupling to gravity that 
appears in the effective action of the twisted theory can be written (6],g] in terms of data 
of the Seiberg-Witten elliptic curve 


X, (Ma\ a /A 

l0g UJ + 8 l0g (A4 


A \ 




(3.3) 


where a{u ) is the a-period of the Seiberg-Witten differential, u the coordinate on the 
moduli space, and A is the discriminant of the curve. This should agree with ( |3.2| ) on a 
compact four-dimensional hyperKahler manifold. There are two such manifolds: T 4 and 
K3. For T 4 , x — a — 0, and both ( fTT2| ) and ( |3.3| ) are trivial. For K3, x — 24, a — —16, 
and by comparing ( |3.2|) to (|3.3|) we find: 




(3.4) 


Using the explicit expressions A = u 2 — A 4 , and the explicit relation between u and a 
worked out in |24|], 


1 


u(a) — 2a ( 1+24 + + FT 2 \TI + 


5 /A\8 9 /A\ 12 1469 /A\is 


2 12 V i 


y (A 
W Y 


+ 


2 2S 


-) + 

a. 


(3.5) 


we computed the first few non-trivial corrections: 


6 


A 


a 


wn, x 1 / A \ 8 1 / A\ 12 1647 / A \ 16 

^)( O )=-l 0g (32 T ) +5S (-) +3^n(-) + 


1647/Ay 


981 /A\ 20 450137 /A\ 24 45111 /A\ 2 8 


(3.6) 


2 42 V 


a 


2 31 V a . 

This is precisely what we found in sect. 2. 

It is interesting to notice that, when JF- 1 ) is written in terms of modular forms (by 
using for example the results in H), one finds 


JF (1 ) — _ log 77(7-). 


(3.7) 


In this equation r is the modular parameter of the Seiberg-Witten curve when written in the 
r(2) description of 0, and it is related to the modular parameter tq of the r°(4) description 
of [g,[24| by r = To/ 2. Equation ( p.7|) is in agreement with the general considerations in 
PI relating F { 11 to //(r) whenever the local Calabi-Yau geometry reduces to a Riemann 
surface. 
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4. Matrix model considerations 


The conjecture of Dijkgraaf and Vafa states that the cubic matrix model captures all 
the information about the of the II geometry. On the other hand, one can consider 
a limit that recovers the Seiberg-Witten exact solution to SU(2). Af — 2 super Yang-Mills 

], which are given by the non¬ 


theory, together with its gravitational corrections |T7| , 
planar sector of the matrix model. In this section we give some results for the cubic matrix 
model beyond the planar approximation. We compute the free energy up to sixth order in 
perturbation theory, providing in this way information up to g — 3, and we use the loop 
equations to find an exact expression for F i - 1 \ 


4-1. Perturbative calculation 


The cubic matrix model has the potential 


W($) = tr(y$ 2 + |$ 3 


N 

£■ 

i —1 


m x2 g n 

+ T< 


(4.1) 


Classically the potential is extremized if the eigenvalues are in either one of the two critical 
points, i.e. at a± — 0 or a 2 = —m/g. Quantum-mechanically the eigenvalues form two 
bands and there is eigenvalue tunneling between the two bands. We want to consider the 
metastable vacuum in which N± eigenvalues are at 0 and IV 2 eigenvalues are at —m/g, with 
Ni and fixed, subject only to the condition that N\ + N 2 — N. This corresponds, in 
matrix model terminology, to a two-cut solution. This issue has already been addressed 
recently |l6j]. Here we will use a slightly different approach which avoids introducing 
ghost degrees of freedom. While the authors of [[L6| considered only the planar limit of 
the matrix model, we also compute non-planar contributions. While the former contain 
the information about the effective superpotential, the holomorphic couplings of the field 
theory to gravity are obtained from the non-planar part of the free energy. 


The partition function Z and free energy F of the matrix model are: 


Z = W = 


Vol{U(N)) 




N\( 2F 


N 


f ndVA 2 (A)e- ! f£< x *-*E i y 

(4.2) 


13 












where A(A) = Y\j <:j (A, — A j) is the Vandermonde determinant. We expand around the 
vacuum with \ = 0 for i = 1, ..., N\ and \ for i = N\ + 1,... , N . Denoting the 

fluctuations by //., and , the Vandermonde determinant becomes 

a2 ( a )= n n ( u h-vh ) 2 n ■ ( 4 - 3 ) 

1<21<Z2<A/"i 1<.7i <j2<N2 

_ _ 1 <3<N 2 


We also expand the potential around this vacuum and get 


N 1 


*r = E 


2 = 1 


9 3 


m r 

2 14 + V 


N 2 

E 

2=1 


(m 2 

V 2 Vi 


9 -uf 
3 1 


m 

+ W 2 ' 


(4.4) 


Notice that the propagator of the fluctuations around —m/g has the ‘wrong’ sign. The 
interaction between the two sets of eigenvalues, which is given by the last factor in ( |4.3| ) , 
can be exponentiated and included in the action, as in [|18|. This generates an interaction 
term between the two eigenvalue bands 


W int = -2N 1 N 2 \og 


m 


9 


w 1 

k=1 


‘EE 

i,j P =0 


(-1) ? 


P 


p k—p 

■ 


(4.5) 


By rewriting the partition function in terms of matrices instead of their eigenvalues, we 
can represent this model as an effective two-matrix model, involving an N± x N\ matrix 
$1, and an N 2 x N 2 matrix $2: 


Vol(f7 (Ni)) x Vol (U(N 2 )) 


T»$iT»$ 2 e _w " l($l)_W2(<i>2)_w($1 ’ $2) , 


(4.6) 


where 

= + tr(im$?+ !$?), 

W 2 (4> 2 )=-tr(im$2-|$3j 5 

00 1 k /i\ (4-7) 

Winter, $2) = 2V-(lf V(~1) P tr tr $ 2 _P 

fci k m U i vw 

+ N 2 W(a 2 ) + N]W (aj.) + 2 N ± N 2 ln(^) ■ 

where tr^ = Ni, tr<f >2 = W 2 , W(ai) = 0 and W(a 2 ) — We have dropped the statis¬ 
tical factor NI/(Ni\N 2 \) which counts the number of ways to distribute the N eigenvalues 
among the two critical points of the potential. This two-matrix model is (perturbatively) 


14 



well defined if we choose <f>i hermitian and $2 anti-hermitian, i.e. ji, real and Uj imag¬ 
inary. We can now compute the free energy F — log (Z) is a straightforward way. F 
consists of two parts, a perturbative part F pcrt = \og(Z(g)/Z(g=0)^ which vanishes for 
the free (Gaussian) model and a non-perturbative part F n _ p _ which gets contributions from 
the U(N) group-volume factor. Both can be computed in a straightforward way. The 
perturbative part can be expanded as 

00 2 h 

F^ = -N 1 W(a 1 )-N 2 W(a*)+2N 1 N 2 ]n(-')+'%2 ^ (£ 3 ) F M (AT,iV 2 ) (4.8) 

^ h =1 9>o m 

h+2-2g>0 


where Ff l fj is a homogeneous polynomial in Ap and N- 2 of degree h + 2 — 2g. One finds, 
up to h — 6 , 


771 

Tpert = -JW(oi) - N 2 W(a 2 ) + 2 N ± N 2 ln( —) 

. - 5 N?N 2 + 5NiN% - InI) + i(W - N 2 ) 

LV3 3/6 

|w 4 - + 59 Af 2 IVf - y ivwl + |w 4 ) + (l^i 2 - y AhAh + |fV 2 2 


4 




.6 r 


+- 




m 


,10 


+ - 


m 


,12 


m 


56 Ar5 871 t4ikt 2636 at3 Ar2 2636 at2a73 871 Ar at4 56 at5 

—AT-AT Ah H-AT Ah-AT Ah H-Ah Ah-Ah 

3 1 3 1 3 12 3 1 2 3 2 3 2 


Y 332 3 _923 


V 9 


Ar 2 A j 923 A7 A7 2 332 3 \ 35 , A _ AT , 

2 Ah Ah 4— —NiN 2 - —N 2 j+—(Ni — N 2 ) 


12 


f 512^6 _ 6823 


LV 3 


N?N 2 + 2 ^N A l N 2 2 - ± • •.) 

2 3/ 


67310 


/1864 4 47083 at3a7 iro ^ AT 2 Ar 2 

+ (———Ni — ———Ah Ah T 15349 Ah N 2 4- .. 


6 


+ 338A/f - 1632WW + 338Ah 


15 


9152 


W 7 - 45118 Ah® Ah + 247980W 5 W 2 - 540378Ah 4 Ah 3 ± • • • 

w 3 - 120880Ah 2 Ah 4 = ...) 


/ 54416 
V 5 
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3 

/65536 
LV 3 


/ GG 1 QO 

Ni - 187528 N?N 2 + 570066Ah 3 Ah 2 ^ ) + f— 


at7a7 13258178 Ar6Ar2 37761034 Ar5 Ar3 52780010 Ar4 Ar4 

N{N 2 +---ATW 2 ---AT Ah 3 +---AT AT =f 


-(Ah ■ 

- Ah) 


Ah 8 - 

1933906 


(1762048 6 12980560 at5 Ar 54863776 Ar4 at2 256344964 3Ar3 , 

V 9 1 3 1 3 9 12 

/1305280 Ar4 18059582 A o „ , „„ Ar2 Ar2 

+ ---AT---AT Ah + 11824166ATAT 4= . 


+ 


71680704 ,, r2 2 \ 8748896 


-(Nt + Nj) - 


JV 1 JV 2 


+ ... 


(4.9) 
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where ... means again anti-symmetric completion. The non-perturbative contribution to 
the free energy is 


e Fn. p. = m -H N ?+ N i) 


(27r) 


~^(n 1 +n 2 ) 


Ni 


N 2 


n r w n r w- 


(4.10) 


k= 1 


1=1 


With the help of the asymptotic expansion (see e.g. 

N . j\t2 i 

\r(k))= -In AT -In N- , 

2 12 4 2 


, ^ 7y 2 14 1 °° 

ln ( II r ( fc )) = — In IV - - In IV - -iV 2 + -iV In 2^ + C'(-l) + 5^ 


k =1 


this becomes 


B 


2a 


Ma - i) n 2 9~ 2 

(4.11) 


F - = 5 JV * ln ® + 2 


ln ((h) _ |(JVf + Nl) - ± ln( Nl l V 2 ) 


+ 2C'(-1) + ^ 

3=2 


B 


2» 


^ 4.9(.9 - 1) V^- 2 


1 1 

+ 


N. 


2g-2 


(4.12) 


Comparing with the prepotential of the gauge theory with cubic superpotential, which can 
be extracted from the results given in |19[], we find it to be in precise agreement with the 
leading terms, at each order in the coupling constant, of the free energy computed from 
the matrix model if we identify iVj —> Si. The only part of the gauge theory result which 
is not directly determined by the matrix model is the dependence on the cut-off A. 


4-2. F W from the loop equations 

The expression for the free energy ( |4.9| ) was obtained by doing perturbation theory. 
However one can obtain analytic expressions for the genus g free energies by using the 
loop equations of the matrix model. These equations were used systematically in [31] in 


the one-cut case, and extended to the multi-cut case in [32], Before moving on to the 


comparison with topological string theory calculations, we derive an analytic expression 
for for the cubic matrix model which can be easily expanded to high powers in Aq 
and N 2 . The computation closely parallels the derivation in |[T^| of the two-cut solution, to 
which we refer for further details. The reason why we cannot directly adopt the result of 
|3^1 is that instead of imposing the absence of eigenvalue tunneling between the two bands 
(equal chemical potential), we fix Aq and N 2 , i.e. we impose 


r*X 4 


p(x)dx — Ni , 


cx 2 


and 


p(x)dx — N 2 . 


(4.13) 


'*3 


' X 1 
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instead of p(x)dx — 0. In other words, we are not looking for the true vacuum but 
consider a meta-stable vacuum with fixed N± and Here p{x) is the eigenvalue density 
which is given by the discontinuity of the resolvent of the matrix ( I>. For the cubic superpo¬ 
tential W{x) = ^x 2 + §x 3 it is p(x) = V / ft^) = ^iV w>2 ( x ) ~ h( x ) where 
/1 is a polynomial of order one whose two coefficients parameterize the widths of the two 
branch cuts. Note that Ni are given by the same integrals as S r in the gauge theory. The 
two conditions ( f4.13|) are not independent since IV 1 + IV 2 = N. The other conditions, which 
follow from the asymptotic behavior of the resolvent, are exactly as in P^fl. Following the 


steps in [|32f| one finds (we give a few intermediate results of the computation in appendix 

B) 


1^11 1 1 

T (1) = y In Mi --In K(k)~— ^ \n{xi-Xj) 2 +- ln(xi-x 3 ) 2 + - ln(x 2 -x 4 ) 2 +const. 

i=l i<j 

(4.14) 

where K(k) is a complete elliptic integral with modulus 


2 (xi - x 2 )(x 3 - x±) 


k = 


(xi - x 3 )(x 2 - x 4 ) 


(4.15) 


and the moments Mj are defined as 


M, = 


dx 


W'{x) 


2 ni 




(4.16) 


where the contour C encloses both cuts. For cubic W, Mi — g for all i — 1,... ,4, as one 
shows by deforming C such as to enclose infinity. 


It remains to express F B) as a series in N± and N 2 , with coefficients depending on m 
and g , the two parameters of W. For this purpose it is convenient to change variables, as 
in |1, (xi,x 2 , x 3 , X 4 ) —> (A 2 i, A 43 , Q , /), where the explicit relations were given in ( |2.18|) . 


4 The complete matrix-model partition function involves a sum over all possible eigenvalue 
distributions and does not have a topological expansion, as explained in [ff3|. This subtlety is, 
however, not relevant here as we fix Ni and N 2 . 
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Inserting ( 2. 19|) with S', replaced by Ay. we obtain 

g G f 332Nf 923N?N 2 , \ g 8 (1864N? 47083Ab 3 A 2 ir0/1 „ Ar2Ar2 

+ f - -± • • • ) + I -^- + 15349Ai 2 ./V 2 2 ip 


V 9 


m 12 \ 


6 


+ -^-nr ( 5441 ^ 6iVl _ 187528 A^i 4 A^ 2 + 570066A^i 3 77 2 2 T • • • 


+ 


m 15 \ 5 

# 12 /1762048Ab 6 12980560 A^i 5 77 2 54863776A^i 4 77 2 2 256344964Ab 3 N 2 3 


TO 18 \ 

+ 0(N 7 ) 


9 


+ 


9 


± ... 


(4.17) 

where . .. means, as before, antisymmetric completion. This expression agrees exactly 
with the expansion of for the local geometry II (|2 .2 1|) that we obtained in section 
2, with the identification S’, — Ni, i = 1,2. This provides a one-loop test of the relation 
conjectured in |[bj . 


Using the iterative procedure developed in [fTy, one can, with some effort, derive 
expressions for the higher genus contributions to the free energy. 


4-3. Af = 2 Yang-Mills from the matrix model 


As explained in [T^] |[Uij| ||TG[ . using the results of the cubic matrix model one can derive 
results for SU( 2), J\f — 2 Yang-Mills theory. The idea is the following: the cubic matrix 
model corresponds to the Af — 2 theory broken down to Af = 1 by adding the tree level 
superpotential IU($) for the chiral superfield ( I>. By taking g. rn —> 0 and keeping A fixed, 
we recover the pure Af — 2 theory. In general, a tree level superpotential of order r + 1 
allows one to recover the moduli space of the SU{r ) theory. In our case, we can recover the 
SU( 2) theory, and A will then be related to the well-known 'u-modulus of Seiberg-Witten 
theory 0H. Some results for the prepotential for general r have been recently obtained 

a- 


m 


In order to recover the Af — 2 results, we first have to find the relation between the 
matrix model variables and the usual variables of the Seiberg-Witten solution. Following 
0 , we extremize the effective super potential of the gauge theory TTy.fr = dSi 


d F (0) | <9F (0) 

dS 2 
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where (cf. (4T) and ( 4.12 ) with iVj —> 5* and the A dependence added on dimensional 
grounds) 


jr(°) 


1 

2 


+ log (^) - i (s2 + s2 > + 23 A 

+ (jS? - 5SfS 2 + 56'i5| - jS|) + 0(S 4 ). 


(4.18) 


Here we have ignored terms 0(S ) since they will not contribute to the positions of the 
extrema of W e ff- Extremizing W e fj, he. solving ds t W c ff = 0, requires Si — —S2 = S with 


— = f-l 4 

g A 3 VA ) 


,A\4 /A\8 / A \ 42 /A\16 / A\20 

1+6( a) + 140 (a) + 4620 (a) + 180180 (x) + 7759752 (^j +••• 

(4.19) 


as obtained in |T6[]. Note that the relation S\ + S 2 = 0 implies that fi in (| 2 . 2 |) vanishes. 
This, on the other hand means that the contours which define Si and S 2 can be deformed 
into each other without picking up a contribution from the point at infinity. In other 
words, we can add this point and consider a compact Riemann surface, as in Seiberg- 
Witten theory. However, the scale A which appears in ( |4.19|) is, a priori, not identical with 
the scale appearing in Seiberg-Witten theory. To find their relation, we relate the curve 
y 2 = W'(x) 2 + fo — 9 2 II(* — x *) fh e Seiberg-Witten curve. First we shift x —> x — ^ A. 
Using the relations (|2.18| ) and the solutions A 2 i(S, —S) and A 43 (S, —S) (c.f. ( |2. 19| ) ) one 


finds, after rescaling y —> gy, y 
SW curve, y 2 — {x 2 — u) 2 


x‘ 


-iA : 


— 4A 4 . Comparing this with the SU( 2 ) 


Ag W leads to the identihcations 


1 A2 

“ = i A ' 


and A = —— Agw 

V2 


(4.20) 


This, together with ( [4.19Q , gives the relation between the variables in the matrix model 
computation and the usual variables in Seiberg-Witten theory. 


Since in order to obtain pure 577(2) super Yang-Mills we take a limit of the matrix 
model, in order to compare the results of the matrix model calculation with the results 
obtained in sections 2 and 3 we have to be very careful and look for quantities that do not 
vanish or diverge as g —> 0. For example, 5 vanishes as g —> 0, while S/g is independent of 
g and can be expressed solely in terms of A, as it is apparent from ( |4 .1 9| ) . When we express 
F( r ) in terms of S/g , we see that it is given by g 2 2r , times a function of S/g. Therefore, 
if r 7 ^ 1, the resulting quantity depends on g and vanishes or diverges as g —» 0. This is 
related to the fact that F ( C is not a function in moduli space, but rather a section of a line 
bundle £ 2-2r , i.e. for r 7 ^ 1 it is not invariant under the gauge transformations of special 


19 






















geometry. This also indicates that, in order to compare the matrix model calculations with 
the results of section 2 and 3, we should focus on gauge-invariant quantities. 


Following this idea, in the case of F ii]] it is clear that the appropriate quantities are 
the second derivatives w.r.t. the moduli, i.e. the r v/ couplings, which are independent of 
One can check 


g |1_£U - 'trie can check [lb| that the T t] computed in the matrix model lead to the right 
result for r in the AT = 2 theory. Since is gauge-invariant and can be expressed in 
terms of A, with no g dependence, one should be able to compare directly the F ^ of ( |3.6p 
with the F l 1 ' obtained in the matrix model. Indeed, using ( |4.19| ). ( |4.20|) , and the relation 
between the u modulus and the a variable given in (|3.5|), we find that the F W of ( [4.17] ) 
reproduces ( |3.6|) after evaluating it at the extremum I. It seems less straightforward to 
relate the matrix model results for F^ for g > 1 to the couplings JFO) j n Seiberg-Witten 
theory derived from string theory. In the final section we will, however, find agreement 
between the expressions ( |2.17|) and a recent conjecture by Nekrasov. 


5. Comparison with instanton computations 


The semiclassical expansion of the J\f — 2 prepotential of Yang-Mills theory can be 
obtained by direct computation of instanton corrections. In a recent tour de force Nekrasov 
m (see also |ll|]) was able to provide general expressions for the n-tli instanton contribution 
to the J\f — 2, SU(N) prepotential, with or without matter. The answer for the n-th 
instanton correction has the form jF n (a, e 1 , 62 ) and it is an analytic function in e i. 2 - For 
Ci = —62 = e, one can expand this as 

OO 

e~ 2 F n (a,e) = ^^ s) (a)e 29 “ 2 . (5.1) 

k =0 


The first coefficient in this expansion, F^\a), gives the prepotential of the corresponding 

(9), 


J\f — 2 theory. It was conjectured in JHJ that the remaining coefficients Fn (a), 9 > 1, 
give the n-instanton correction to the gravitational couplings of the J\f — 2 theory. 
Here we provide nontrivial evidence for this conjecture by comparing some of the instanton 
computations of HJ with the results of section 2 (and, therefore, with the matrix model 
computations). 


Note that the scale appearing in sect. 3 is, in fact, Asw- 
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In [llj] T n (e 1 , 62 ) was explicitly computed for SU( 2 ), up to n — 4. To compare with 
our results in section 2, we restore units as follows: T n —> (A/2) 4 n jF n . For n — 1 one finds 


e- 2 ^i(e) = 


,-2 


A 4 


2 5 a 2 ' 


(5.2) 


This says that the one-instanton contribution to vanishes for g > 0, which agrees with 
( P.16| ) and ( |2.17|) . For n — 2, one finds: 


e~ 2 Me) 


A 8 10a 2 — e 2 
2 8 8e 2 a 4 (4a 2 — e 2 )" 
5 A 8 I A 8 


_ _ 2 o iy ( 1 iv 2 11 A 8 

+ ^ 13^8 +£ 218 ^10 


+ e 4 


7 A 8 

2^ ^F 2 


(5.3) 


+ 


For n — 3, one finds from |lT 

A 12 


e~ 2 ^{e) 


18a 4 — I3a 2 e 2 + e 


2 12 24e 2 a 6 (a 2 — e 2 ) 2 (4a 2 — e 2 ) 2 


-2 


3 A 12 


I A 12 2 117 A 12 4 

+ x xn To" + + e 4 


2 18 a 10 ' 3 • 2 14 a 12 


2 22 a 14 


293 A 12 

^23 ^16 


(5.4) 


H- 


Finally, for n — 4, one can extract from |1T 
e _2 JF4(6 


A 16 23504a 10 - 70872a s e 2 + 6746fa e e 4 - 26339a 4 e 6 + 3708a 2 e 8 - 162e 10 


2 16 


128e 2 a 8 (4a 2 - 9e 2 ) 2 (a 2 - e 2 ) 2 (4a 2 - e 2 ) 4 


-2 


1469 A 16 1647 A 16 

+ 


2 31 a 14 2 29 a 16 


+ e 


,171201 A 16 
2 34 ^ 


+ e 


985823 A 16 
^3^ ^20 


+ ... 


(5.5) 

We see that, after relabelling e —> it. the coefficients in the above expansions are in perfect 
agreement with the results in (|2.I4|), (|2.16|) and ([2.1 7|) . 
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Appendix A. 


We collect the first few polynomials gl; l> which enter the topological amplitudes (| 2 . 12 |). 


hf ] = (1 - «) 4 

h'.g = 1(2 - nq + 80 q 1 + II 9 O 5 3 + 8 O 5 4 - 17/ + 2/) 

hf ] = - (3 - 38 q + 1712 q 2 + 48326 q 3 + 124634 q 4 + ...) 

3 

hf } = i(4 - 67 q + 14592 g 2 + 711059 q 3 + 4924138 g 4 + 9244668 q 5 + ...) 

hf ] = - (5 - 104 q + 76705 g 2 + 6090098 q 3 + 82568187 q 4 + 358139062 g 5 + 580752958 q 6 + ...) 

5 (A.l) 


d 3) = (i - «) 6 

h< 3) = 1(1 + 4g + /)( 8 - 125j + 1016/ - 8854/ + 1016/ - 1255 8 - 85 ®) 

8 

h {3) = -1(3 3 - 428 q + 207 q 2 - 584608q 3 - 6606954g 4 - 13969512g 5 - ...) 

O 

h i3) = ^ (4 3 - 1275g - 62364g 2 - 14760656g 3 - 289184988 q 4 - 1451906781g 5 - 2450425504g 6 - . 

h {3) = ^(5 3 - 2994g - 655802g 2 - 183946424g 3 - 5822909304 q 4 - 53583541710<? 5 
o 

- 190175648587g 6 - 288031344528g 7 - ...) 

(A. 2 ) 

They were used to derive ( |2.17j ). Note that the can be used to write down the 
conjecturally integer Gopakumar-Vafa invariants for 0(— 2 , — 2 ) —> F 0 (see e.g. prf ) n^ d2 
for g < 4, d\ < 6 and d ,2 arbitrary, e.g. 

{n { i\,d 2 = 0,...} ={0, 0,15, 4680,184056, 3288688, 36882969, 300668486, 

1935031484,10359890196, 47820549652,195274337280, 719145083800,...} 

We checked integrality of the n^ d2 for g < 4, d\ < 6 and d 2 < 2000. 


Appendix B. 


We collect some of the intermediate results which are needed to derive (|4.14|) . As we 
mentioned in sect. 4, the derivation is almost identical to that in sect. 6 of Jj2|. Instead of 


22 








the elliptic integrals Ki defined there, we encounter (we choose x 4 < x 2 < x 3 < x 4 ) 


K, = 


rX4 y/(x — Xi)(x — X2)(x — X 3 )(x — X 4 ) 


'x 3 


(X ~ Xi) 


dx 


(B.l) 


They arise from varying the constraint that N\ = j X4 p(x)dx is kept fixed rather than 
requiring JJ ' 3 p{x)dx — 0 as in [j32j . Explicit calculation gives 


Ki — zX{(—a 6 + 2cr 4 + a 2 k 2 — 2 a 4 k 2 )E(k) — (k 2 — a 2 ) (a 4 — 2a 2 + k 2 )K(k ) 

+ (a 8 — 4cr 6 /c 2 + 6 a 4 k 2 — Aa 2 k 2 + A; 4 )n(a 2 , k 2 ) j 
K 2 = iX{-a 2 (a 4 - 2cr 2 - 2 a 2 k 2 + 3 k 2 )E{k) - ( k 2 - a 2 )(a 4 - 2a 2 + 4a 2 A; 2 - 3 k 2 )K{k) 
+ (a 8 — 6a 4 k 2 + 4a 2 k 4 + Aa 2 k 2 — 3A; 4 )n(a 2 , A; 2 )} 

K 3 = iX{(3a 6 - 2a 4 - 2 a 4 k 2 + a 2 k 2 )E(k ) + ( k 2 - a 2 )(3a 4 - 2a 2 - k 2 )K(k) 

+ (—3a 8 + 4a 6 + 4a 6 A; 2 — 6a 4 A; 2 + A; 4 )n(a 2 , A:)} 

K 4 = iX{(-a 6 - 2a 4 + 2a 4 k 2 + a 2 k 2 )E(k ) - {k 2 - a 2 )(a 4 + 2a 2 - 4a 2 A; 2 + k 2 )K(k) 

+ (a 8 — 4a 6 + 6a 4 A; 2 — 4a 2 A; 4 + fc 4 )II(a 2 , A:)} 

(B.2) 


where 

2 _ (x 4 ~X 3 ) 

(x 4 - X2) 

k 2 _ fol ~ X 2 )(x 3 ~ X4) 

(Xi - x 3 )(x 2 ~ X4) 
x _ 1 (xi - x 3 ) 3/2 (x 2 - x 4 ) 7 / 2 
4 (x 4 - x 3 ) 2 (x 3 - x 2 ) 


(B.3) 


Next we have to solve 
above. We find 


eqs. 


( 6 . 1 ) in 




for the coefficients a^, however with K r as given 


a 1 


a 2 


«3 


a 4 


1 


(xi 

- x 2 ) 


1 

(x 2 

- Xi) 


1 

(x 3 

- X4) 


1 

(x 4 

- X 3 ) 


(x 4 - x 2 ) E(k ) 
(#4 — xi) K(k) 
(x 3 - X\) E(k) 
(x 3 - x 2 ) K(k ) 
(x 4 - x 2 ) E(k) 
(x 3 - x 2 ) K(k) 
(x 3 - xi) E(k) 
(x 4 — xi) K(k) 


(B.4) 


With these values for k 2 and a,, the derivation of F ( 1! proceeds exactly as in ||32|| and 
finally leads to (|4.14|). 
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